Processing Analog Signals With Respect to 
Timing of Related Rotary Parts 
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Priority Data 
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serial number 10/039,193, filed December 31, 2001, entitled 
"Vibration Tracking Filter", which claims the benefit of U.S. 
provisional patent application serial number 60/259,353, filed 
December 30, 2000, entitled "Vibration Tracking Filter". 

Technical Field 

This invention relates to processing analog signals such as 
vibration-indicating signals by filtering, integrating, converting to 
digital signals in response to said rotary signals and resolving into 
orthogonal components, of which the sum of the squares are 
determined and averaged. Pulse frequency multiplication, shared A 
to D, start of averaging, determining sines and cosines of non-integer 
pulses per revolution, and built-in test using square waves of which 
the higher harmonics are filtered out are involved in the processing. 

Background Art 

When an analog signal is to be sampled for A/D conversion in 
response to signals received from a rotating device, so as to permit 
resolving the analog signal into its orthogonal components with 
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respect to the instantaneous angle of the rotating device, aliasing 
may occur near harmonics of the sample frequency. In order to avoid 
that, the timing sampling signals may have the frequency thereof 
multiplied so as to avoid the aliasing. However, frequency 
5 multiplication circuits that can multiply signals of a randomly varying 
frequency and different phases, for example, phase lock loops, are 
complex and require extensive circuitry to implement. 

When processing signals through a multiplexer, signal 
conditioners and an A/D converter, it is most economical if a single 

10 channel can be used to process signals at a plurality of different 

timing sequences. If the timing sequences are at mutually unrelated 
frequencies and have no set phase relationship, many A/D 
conversions may be made without conflict between them. However, 
when an A/D conversion has been started by one stream of control 

1 5 signals, and another A/D conversion is requested by another stream 
of control signals, a conflict exists; if the conflict is not resolved, a 
set of data can be lost. 

When signal conditioning circuitry, including bandpass 
filtering and integrating, feeding an analog to digital converter, are 

20 shared by a plurality of input signals, which are selected by a 

multiplexer, there is a need for a settling time when the multiplexer 
switches from one input to another input before the output of the 
signal conditioning circuit fairly represents the selected input signal 
due to the need to allow the integrator and bandpass filters to settle 

25 on the new value. If the sampling rate of the A/D converter is not 
constant, and particularly, if it may become slower as the result of 
deceleration of a moving part or other factor, the amount of time for 
settling, within a fixed processing cycle, may vary. Therefore, the 
amount of time which can be allocated to averaging of signals within 
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that fixed cycle may vary in dependence upon whether the timing 
sampling signals are slowing down or not. 

To resolve components of a signal into orthogonal 
components at angles of a related rotating part, it is common to 
5 utilize lookup tables of the desired functions in response to the angles 
of the rotating part. However, if the rotating part has a non-integer 
number of pulses per revolution, table look-up does not seem to be 
appropriate. 

When signals which are processed from sensors of real time 
10 phenomena are used either to control or to monitor the safety of 

equipment in use, and particularly equipment in use in flight, such as 
aircraft jet engines, correct operation of the signal processing 
apparatus must be assured. For this purpose, it is known to utilize 
built-in test equipment (BIT) to periodically test critical signal 
1 5 processing channels to determine the viability thereof. When the 
processing involves the frequency and amplitude of AC signals, the 
utilization of sine waves as test signals may be involved. Sine waves 
may be generated by oscillators, of course, which requires a fair 
amount of circuitry per oscillator, most oscillators being operable only 
20 at one frequency, or possibly at several frequencies with auxiliary 
apparatus. When many frequencies are involved, many oscillators 
would be required. Sine waves may be synthesized stepwise and 
smoothed, but this requires significant digital processing, particularly 
for a large number of frequencies that may be used. 



25 Disclosure of Invention 

Objects of the invention include: providing simple 
multiplication of pulse frequencies on a rapid real time basis; allowing 
signals sampled at two different rates to utilize the same filtering, 
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integrating and A/D conversion channel, without loss of samples due 
to conflicts between the diverse sampling rates; determining the 
maximum time which is required for signal processing, thereby 
allocating as much time as possible to settling of signals, after the 
5 input to bandpass filter and integration circuitry has changed from 
one input signal to another; provision of a simple table lookup of 
trigonometric values of angles of a rotating part which provides a 
non-integer number of pulses per revolution; and provision of a 
simplified source of test signals for testing sinusoidal signal 

10 processing channels involving filters. 

This invention is predicated in part on the realization that 
pulses may be inserted between the pulses of a stream, so as to 
multiply the pulse frequency thereof, utilizing the period between a 
given pulse and the next preceding pulse as a determination of the 

15 length of each fraction of a period when a pulse should be inserted. 
This invention is further predicated on the concept that conflicting 
demands for A/D conversion of filtered and integrated analog signals 
need not result in loss of samples if the conversion for one of the 
timing signals is used as the conversion for the conflicting timing 

20 signal. The invention is predicated in part on the discovery that a 
lookup table of trigonometric values, corresponding to a non-integer 
number of pulses per revolution from a sensor of a related rotating 
part, can be approximated by selecting a different number of pulses 
per revolution which is close to the actual non-integer number, and 

25 using a count of a small integer multiple of that different number of 
pulses to provide repetition of the trigonometric values desired. This 
invention is also predicated in part on the realization that the effect of 
upper harmonics of lower frequency components of a square wave, 
whose magnitude is adjusted upwardly by 4/n to produce a sinusoidal 
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fundamental frequency at the amplitude and frequency of a sine 
wave of interest, is predictable. 

According to the present invention, analog signals, such as 
from a vibration sensor, are processed first through a single channel 
5 of filtering, integration and conversion to digital signals, and then 

through a plurality of channels in which the signals are resolved into 
orthogonal components in accordance with the instantaneous angles 
of related rotary parts having sensors which emit pulses per 
revolution; the pulse frequency is multiplied by inserting bits between 

10 the bits from the sensors; maximum settling time for filters and 

integrators following switching from other signals, such as square 
wave test signals, or other sources, is achieved by starting the 
processing in each separate channel at a point in time which allows 
just sufficient time to accommodate a maximum deceleration in that 

1 5 channel from current speed; if one channel seeks to start the A/D 

while it is processing for the other channel, the digital output is used 
in both channels; non-integer number of pulses per revolution of a 
sensor is accommodated by counting modulo-one-bit-higher than a 
small integer multiple of another number close to the non-integer 

20 number of interest; built-in test passes square waves through all of 
the processing channels in an interleaved fashion with input signals, 
the predictable response of the processing channels to the square 
wave being indicative of either proper functioning or fault in the 
processing channels. 

25 According to the present invention, the period between one 

pulse and the next preceding pulse is divided into N, and N minus 1 
additional pulses are inserted between each pulse and the next 
succeeding pulse, thereby to provide pulse frequency multiplication. 
More than one rate of sampling of an analog input is made by an 
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analog to digital converter; according to the invention, whenever 
there is a conflict resulting from one sampling signal attempting to 
start to an A/D conversion while an A/D conversion is in process, the 
resulting digital output will be utilized not only for the signal that 
5 started the conversion, but also for the signal which is in conflict 

therewith. In further accord with the invention, the portion of a fixed 
cycle required for processing signals, such as averaging a large 
number of signals, which are sampled in response to timing signals 
from a rotary part related thereto, the time required for such 

1 0 processing in the event that the rotary part decelerates at a maximum 
predetermined deceleration rate is determined, and the processing is 
started at a corresponding time advanced from the end of the fixed 
cycle. In still further accord with the invention, provision of 
trigonometric functions of an input signal in dependence upon the 

1 5 angle of a rotary device which provides a non-integer of numbers per 
revolution, is determined by selecting a different number, integer or 
non-integer, of pulses per revolution, counting modulo a number 
which is a small integer times the selected number of pulses per 
revolution, and providing a lookup table of the desired trigonometric 

20 functions responsive to the count. In accordance with the invention 
still further, built-in test equipment which tests the channels for 
processing an analog signal utilizes square waves at an appropriate 
frequency, the third and higher harmonics of which are suppressed 
by the anti-aliasing filters of the signal processing apparatus. 

25 Other objects, features and advantages of the present 

invention will become more apparent in the light of the following 
detailed description of exemplary embodiments thereof, as illustrated 
in the accompanying drawing. 



Brief Description of the Drawings 

Fig. 1 is an illustrative schematic block diagram of a system 
incorporating the invention. 

Fig. 2 is a block diagram of signal conditioning circuits. 

Fig. 3 is a schematic block diagram of a 2 n pulse frequency 
multiplier according to the invention. 

Fig. 4 is a logic flow diagram of a tracking filter state 
machine. 

Fig. 5 is a partial table of sine and cosine values of fan 
angular rotation. 

Fig. 6 is a partial table of sine and cosine values of gas 
generator angular rotation. 

Fig. 7 is a time chart illustrating the determination of 
averaging start times according to the invention. 

Figs. 8 and 9 illustrate a three-pole filter. 

Mode(s) for Carrying Out the Invention 

Referring to Fig. 1 , an exemplary task which illustrates and 
utilizes the various aspects of the present invention is sensing 
vibration in a device such as a gas turbine engine which has two 
engine portions rotating at different speeds, such as a gas turbine 
engine fan and a gas turbine engine gas generator. The signals from 
the vibration sensor are analyzed to determine the magnitude of 
vibration related to each of the rotary speeds. However, it should be 
understood that aspects of the invention can be utilized when 
analyzing the output of a vibration sensor with respect to a single 
rotary frequency. Similarly, the invention may be utilized in a time 
division multiplex fashion to process the signals from two, three or 
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more vibration sensors, or to analyze the output of the vibration 
sensors with respect to more than two rotary frequencies. 

In Fig. 1, a vibration sensor 20 may comprise a piezoelectric 
accelerometer which provides signals to a charge amplifier 21 . The 
5 output of the charge amplifier 21 may be selected by a multiplexer 
22 so that the signal will be processed, or the multiplexer 22 may 
select signals from a square wave test signal source 23, as is 
described more fully hereinafter. 

Whichever signal is selected by the multiplexer is provided 

1 0 over a line 27 to a circuit 26 of filters and an integrator which is 

illustrated in Fig. 2. In Fig. 2, a pair of high pass filters 30, 31 each 
have a cutoff frequency, fc, of 18.3 Hertz, the filter 30 having a 
damping factor, 5, of 0.6 and the filter 31 having a damping factor of 
0.707. A low pass filter 32, which may be thought of as a pseudo 

1 5 integrator, converts the acceleration signal that has been passed 

through the high pass filters into a velocity signal. The output of the 
integrator 32 is passed through three low pass filters each having a 
cutoff frequency of 630 Hertz, with various damping factors as 
shown in Fig. 2. Then an amplifier 36 provides a gain of 61 .1 . 

20 Finally, a high pass filter 37 with a cutoff frequency of 18.3 Hertz 
and a damping factor of 1 .0 completes the signal conditioning. The 
net effect of the low pass and high pass filters is band pass filtering 
of the input signal to avoid aliasing (the low pass filters) and to avoid 
low frequency noise (the high pass filters). The output of the signal 

25 conditioning circuitry 26 on a line 40 is applied to an analog to digital 
converter, A/D 41 (Fig. 1). 

The A/D converter is caused to start in response to signals 
which are multiples of the speed sensor pulses, in this example, from 
the engine fan, Nf, and from the engine gas generator, Ng. In order 
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to avoid aliasing errors, the output signal of an Nf speed sensor 43 is 
applied over a line 44 to a frequency multiplier circuit 45, which is 
illustrated in Fig. 3. 

Referring to Fig. 3, a pulse frequency multiplier according to 
5 the present invention responds to the occurrence of each pulse, such 
as the pulses Nf indicative of fan speed in a gas turbine engine, by 
inserting a requisite number of pulses between one received input 
pulse and the next received input pulse. In the 4x example in Fig. 3, 
a pulse is inserted at one-quarter of the period after a received input 

1 0 pulse, a pulse is inserted at one-half of the period after a received 
input pulse, and a pulse is inserted at three-quarters of the period 
after a received input pulse. The measure of the period, and 
therefore the fraction of a period where a pulse should be inserted, is, 
however, based upon the actual period following the next preceding 

1 5 pulse, which is very nearly of the same period as the present pulse. 

The period of the preceding pulse is manifested in the number 
of counts of a one megahertz clock signal which appears on a line 
50, and is allowed to be counted in a 12 bit presettable counter 51 
only when an AND circuit 52 is enabled by a latch 55. The latch 55 

20 is set by a signal on a line 56 just after a leading edge detector 57 
senses the presence of an input pulse on the line 44. When the 
leading edge detector senses an input pulse, it provides a data 
transfer signal on a line 61 that is fed to a delay circuit 62 which in 
turn has its output fed to a delay circuit 63. The delays in the delay 

25 circuits 62, 63 are on the order of a few microseconds and are 

present just to avoid glitches in the process of setting and clearing 
the latch 55. The data transfer signal on the line 61 resets the latch, 
thereby ending the passage of clock signals from the line 50 through 
the AND circuit 52 into the counter 51 ; that is how the counter is 
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stopped to determine the extent of the period between the present 
input pulse and the next prior input pulse. 

At the same time, the content of the 12 bit presettable 
counter 51, except for the two lowest order bits, is transferred into a 
5 10 bit holding register 66 by the data transfer signal on the line 61 . 
Just after the data is transferred from the 12 bit presettable counter 
51 into the ten bit holding register 66, the output of the delay circuit 
63 will again set the latch 55 so that counting of the present period 
can occur in the 1 2 bit presettable counter. 

10 The data in the counter 66 is truncated by eliminating the 

first and second least significant bits. This is equivalent to a shift to 
the right (toward lower significance) and is equal to division by four. 
Therefore, the third through twelfth bits of the twelve bit resettable 
counter 51 represent the value of one-quarter of the period of the 

1 5 next preceding pulse-to-pulse interval. 

When the latch 55 is set, clock pulses which pass through 
the AND circuit 52 are applied not only to the twelve bit presettable 
counter 51 but also through another AND circuit 67 to cause a ten 
bit up counter 68 to count clock signals derived from the line 50. 

20 The output of the 10 bit up counter 68 is applied to a ten bit 

compare 71, and when it matches the value set in the 10 bit holding 
register 66, a match signal on a line 72 indicates that one-quarter of 
a period has elapsed. This signal passes through an OR circuit 73 to 
become part of the stream of 4x input pulse signals on a line 74. 

25 The signal on line 72 is also provided to a delay circuit 77, the output 
of which will pass through an OR circuit 78 to clear the ten bit up 
counter 68 to make it ready to count the next quarter of a cycle. 
The OR circuit 78 also responds to the data transfer signal on the line 
61, indicating that an input pulse has just been received. 
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The signal on the line 72 is also provided to the clock input of 
a two bit up counter 79, the output of which is applied to a decoder 
81 which will decode binary four (1 1 ) to provide a signal on a line 82 
that will inhibit the AND circuit 67, thus preventing clock signals 
5 from occurring during the fourth quarter of the period. The reason 
being that the fourth pulse will be that provided as the data transfer 
signal by the leading edge detector 57 when the next input pulse is 
received on the line 44. The match signal on line 72 will be provided 
at the end of the first, second and third quarters, and the actual input 

10 pulse, provided by the leading edge detector 57 onto the line 61 , will 
provide the fourth pulse. Thus, there are four times as many pulses 
on the line 74 as are received on the line 44, effectively multiplying 
the input pulse frequency by four. 

In the 4x example herein, n = 2, and 1/2 n is equal to one 

1 5 quarter; 2 n = four. 

The embodiment described in Fig. 3 multiplies the frequency 
of input pulses by four. This is because the data in the 1 2 bit 
presettable counter 51 is truncated by four, by dropping the two 
lowest order bits, meaning that the data is divided by four; this 

20 results in the match signal on line 72 being generated every quarter 
of a cycle. If the output of the 1 2 bit resettable counter was 
truncated by three bits, equal to a division by eight, then the match 
signal on line 72 would be present every eighth of a cycle. Similarly, 
any power of two may be selected, and adjustment of the truncation 

25 of data at the output of the twelve bit resettable counter 51 will 

accomplish the desired frequency multiplication of the pulse stream. 
Thus, multiplying the frequency of input pulses by 2n (which is 4 in 
Fig. 3) is accomplished by truncating to eliminate the n lowest 
ordered bits (2 in Fig. 3). 



- 12 - 

The frequency of a signal on a line 89 from the Ng speed 
sensor 90 is doubled in a frequency multiplication circuit 92, which 
may be the circuit of fig. 3, truncated by only one bit. 

As an example, in this embodiment, there is only one 
5 vibration sensor, but a useful case would have two vibration sensors. 
For instance, one at the front of the engine and one at the rear of the 
engine, both mounted on stationary engine portions. Only one of the 
vibration signals would be processed through the signal conditioning 
filters and the A/D converter at one time. The vibration signal is 

10 processed to see what vibration components there are at the speed 
of the gas generator and to see what vibration components there are 
at the speed of the fan. These speeds are not synchronous or 
bearing any phase relationship to each other, and the amount by 
which the gas generator speed may exceed the fan speed can vary at 

15 all times. As an example herein, we may consider a situation where 
the gas generator speed is on the order of three times greater than 
the fan speed, with no known phase relationship, and with the 
relationship between the two changing continuously. According to 
the invention, instead of having a full processing channel for Nf and a 

20 full processing channel for Ng, a single channel of processing can be 
initiated by either Nf or Ng because there are many times when there 
will be no interference between one pulse and the other. When there 
is a conflict between one A/D conversion being in process when 
another A/D conversion would be requested (but for the conflict), the 

25 problem is solved simply by utilizing the A/D output as the data value 
for the train of signals at one frequency (e.g, Ng) as well as for the 
train of signals at the other frequency (e.g., Nf). This poses no 
problem since the A/D conversion data is the same for either pulse 
rate except for a slight difference, which might occur due to the 
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slight difference in timing when the A/D conversion is requested. It 
should be noted that the sampling of the vibration sensor may occur 
in response to a pulse from either of the speed sensors. By 
definition, if there is a conflict (that both of these sensors are 
5 providing an A to D request at times which are near one another), the 
data at one of those times is going to be very closely the same as the 
data at another of one of those times. The conflict is resolved simply 
by allowing the signal that will load the A/D into the Ng register to 
also load the A/D into the Nf register (or vice versa). In addition to 

10 loading the A/D into both registers 123, 124, the system will 

increment the sequence counters 81 , 95 for both Ng and Nf. The 
A/D complete signal comes a little sooner with respect to the data 
that might have been generated by the second signal. 

In Fig. 1 , the signals Nf x 4 on line 74 and Ng x 2 on line 93 

1 5 are each applied to a latch 110, 111. The output of both latches is 
fed to an OR circuit 1 1 3 which provides a start A/D signal on a line 
1 14. The A/D converter will then begin the process of converting 
whatever analog value was present at its input, into a digital value in 
a known fashion. Either latch 110, 111, will remain set until a signal 

20 from a corresponding AND circuit 116, 117 responds to a 

concurrence of the set condition of the latch and an A/D complete 
signal on a line 1 20 to provide a related load A/D to Nf register signal 
on a line 121 or a load A/D to Ng register signal on a line 122. Each 
of these signals will reset the corresponding latch 110, 111. 

25 If one of the latches, for example, 1 10, has been set and the 

A/D converter has been started, the A/D converter will be in the 
process of converting a signal related to Nf . If before it finishes, a 
signal is received on the line 89 so that the latch 111 becomes set, 
another start A/D signal will appear on the line 1 14. However, this 
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will elicit no response in the A/D converter. The latch 1 10, however, 
will remain set. Thus, when the A/D completes conversion of the 
signal relating to Nf, the A/D complete signal on the line 120 will 
pass through both AND circuits 116, 117 and cause the digital value 
5 of the vibration sensor signal converted in response to the Nf signal 
to be loaded into both the Nf register 1 23 and the Ng register 1 24. 
In short, when the timing of the fan speed and gas generator speed 
conflict by being within an A/D conversion time of each other, the 
first signal causes an A/D conversion, the digital signal will appear in 

10 the stream of digital signals relating to the second signal, but at a 

timing determined by the first signal, which will be advanced in time 
by a small amount from one that would have been converted in 
response to the second signal. It has been shown that such slight 
advance in a conflicting signal has no effect on system performance. 

1 5 This aspect of the invention will work with signals which are 

of the same or different, fixed or variable frequencies, and of fixed 
phase relationship or having no phase relation, so that conflicts can 
exist; however, it is most valuable when used with a pair of signals 
each of which can vary in frequency and phase independently of the 

20 other. 

A tracking filter 130 is shown in analog fashion within the 
dash lines in Fig. 1, and is illustrated in a manner in which it may be 
performed digitally, in an application specific integrated circuit, by 
means of the state machine shown in logic flow format in Fig. 4. In 
25 Fig. 4, the routine may begin at a start entry point 133, and a first 
test 1 34 determines if a start process signal has been received, 
indicating that the processes should begin. Until the start process 
signal is received, the routine cycles on the test 1 34. 

When the start process signal appears, an affirmative result 



- 15 - 

of test 1 34 will reach a test 1 35 to determine if a new value of 
vibration signal has been entered into the Ng vibration register 1 24 
(Fig. 1) or not. This is so when the latch 1 1 1 is set. If so, an 
affirmative result of test 1 35 will reach a test 1 37 to determine if 
5 either 256 iterations for regular signal processing, or 4 iterations for 
built-in test of the Ng channel, have been completed. If not, a 
negative result of test 137 reaches a subroutine 140 to multiply the 
content of the vibration register for Ng times the sine of Ng. Then a 
three pole filter subroutine 141 is performed on the product, in a 

10 manner described more fully hereinafter with respect to Figs. 8 and 
9. Then the filtered product is squared in a subroutine 142. Similar 
subroutines 145-147 provide the square of the filtered product of the 
vibration register Ng times the cosine of Ng. A subroutine 149 then 
provides the sum of the squares. 

15 Reference to the lower central portion of Fig. 1 illustrates in 

analog fashion, with the same reference numerals, the functions 
which have just been described. 

In Fig. 4, a test 151 determines if there is a signal, described 
hereinafter with respect to Fig. 7, indicating that it is now the proper 

20 time to begin averaging of the sums of squares derived as just 

described, over a large number of samplings of the vibration sensor, 
for data sampled at a rate related to the speed of the gas generator 
(Ng). If averaging is to begin, then a step 152 adds the recently 
derived sum of the squares into the content of an accumulator and a 

25 sum counter for Ng is incremented in a step 1 53. 

A test 1 57 determines if the apparatus is currently in a test 
mode, or not. This determination is made by cycle control software 
in a computer 1 58 (Fig. 1) which also determines the square wave 
test pattern from the source 23 and the selection to be made by the 
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multiplexer 22, as described hereinbefore. If the apparatus is in a 
test mode portion of a major cycle, a counting number, N, is set 
equal to four in a step 1 60. But if the test mode is not in effect, a 
negative result of test 1 57 reaches a step 1 61 to set N equal to 256. 
5 Then a test 1 62 determines if the sum counter has reached the count 
of N, or not. If it has, that means 256 real samples or four test 
samples have been accumulated in step 152. A step 164 sets an Ng 
done flag, and a step 1 65 divides the amount achieved in the 
accumulator step 152 by 256; this is achieved by storing only the 20 

10 most significant bits of a 28 bit word in the transit buffer. 

Truncating the eight lowest ordered bits is equivalent to division by 
256. Then a step 167 will reset the Ng sum counter, which had 
been incremented in step 1 53, and the routine will revert to a test 
1 69 to see if the latch 1 10 (Fig. 1) is set, indicating a new vibration 

1 5 sample for Nf. 

If a test 1 62 is negative, meaning that the averaging process 
is continuing, a test 1 77 determines if a communication 
synchronization signal has been received from the computer or not. 
If it has, this means that the subcycle of processing Ng has ended 

20 (which will not usually be the case) so that the sum counter must be 
reset and the routine ended. However, during the averaging process 
until the sum counter reaches N, if there is no synch signal, then 
processing of vibration values at the Nf rate can occur. 

If either test 1 35 is negative (meaning there is no new word 

25 to process relative to Ng), or if test 137 is affirmative (indicating that 
processing for Ng is done), or if test 1 51 is negative (indicating that 
the time to begin averaging of the Ng results has not yet arrived), or 
test 1 77 is negative (indicating that a sync termination has not 
occurred, or if the sum counter is reset (in step 1 67), the test 1 69 is 
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reached at the beginning of the processing for vibration values at the 
Nf rate. 

If there is no new vibration value in the Nf vibration register 
1 23, a negative result of test 1 69 reverts the program back to 
5 processing for Ng. But if there is a new value, then a test 1 79 

determines if the Nf done flag has been set yet or not. If not, then a 
series of subroutines 180-189, identical to the subroutines 140-149 
but relating to Nf instead of Ng, are performed. Then a test 191 
determines if a signal has been received indicating that it is time to 

1 0 start averaging for values related to Nf . If it is not yet time, then 

processing reverts to test 1 35 for processing vibration values relating 
to Ng. If averaging can begin, then a series of steps and tests 192- 
217 perform the same functions as steps and tests 1 52-1 77, 
described hereinbefore, except with respect to Nf instead of with 

15 respect to Ng. Then a test 218 determines if Ng processing is done; 
if not, the routine reverts to test 1 35. But if so, the routine will 
terminate through an end point 220. The averaging of Nf values is 
thus interleaved with the averaging of Ng values, and real values are 
interleaved with test values. 

20 Each time the process ends through the point 220, it lays 

dormant until the next time that a start process signal causes an 
affirmative result of test 1 34, which is the beginning of a macro 
cycle in which these tests are performed. The test mode which is 
determined in tests 1 57 and 1 97 can simply be the second subcycle 

25 of a major cycle, the first subcycle being allocated to processing for 
Ng and Nf, respectively. If more speed sensors are used, more 
subcycles will be needed. 

To provide the magnitude of vibration at the various angles 
of the gas generator as it rotates, or at the various angles of the fan 
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as it rotates, the output of the A/D converter is resolved into 
orthogonal components relative to the sine and cosine of the 
rotational angle indicated by the respective speed sensors Nf, Ng. 
These are then filtered to eliminate sum frequencies above 10 Hertz, 
5 leaving only difference frequencies. These are then squared and 
summed so as to provide to the accumulator 1 52, 1 92 the sum of 
the squares of the resolved components, which equals the square of 
one-half of the instantaneous vector magnitude of the vibration 
signal. 

1 0 Because accelerometers are subject to wide ranges of noise, 

which can render any particular reading non-representative of the 
actual vibration, averaging over a large number of samples may be 
desired. However, with the band pass filtering and integration of the 
signal conditioning circuits 26 of Fig. 2, maximum settling time 

1 5 should be allowed after switching inputs to the A/D converter 41 
from test signals to vibration signals, or from one test signal to 
another. The amount of time which is available to provide averaging 
and the amount of time which is available for settling of the signals in 
the time shared signal conditioning of Fig. 2 must therefore be 

20 managed. 

An exemplary situation is briefly illustrated in Fig. 7. When 
the multiplexer switches from a square wave test input back to the 
vibration sensor input, or vice versa, the signal conditioning filters 
and integrator 26 of Fig. 2 will take some settling time to actually 

25 settle at a present value of the vibration sensor or test signal. It is 
only after this adjustment to the vibration sensor or test signal that 
the A/D converter could provide actual digital representations of the 
vibration sensor or test magnitude. Since the definitive settling time 
can normally be established only by making it extremely long, 
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instead, one aspect of the invention is allowing whatever settling 
time there can be while still permitting sampling 256 values of the 
vibration sensor within the allotted time of a cycle, and averaging as 
described with respect to Figs. 1 and 4. If the rotational speed 
5 signals were constant, then a fixed amount of time could be reserved 
for averaging. Since, however, these speeds are not constant, care 
must be taken to allow sufficient time, even when there is a 
deceleration of either of the rotating members which the speed 
sensors Nf, Ng respond to. 

10 The basic processing time for the minimum speed of the fan 

and for the maximum speed of the gas generator are given, only as 
examples, for use in establishing values on Fig. 7. A deceleration 
(decel) of 40% per second will be 13% in 0.32 seconds, the 
exemplary cycle time herein. Worst case allowance is therefore 

1 5 taken to be an increase of 1 3% in required time to average 256 

samples, over the time that would be required at the current speed. 



Nf: 

min speed = 3,420 RPM; 57 RPS 
20 x 1 2 ppr = 684 Hz 

x 4 (Fig. 3) = 2736 Hz 

period per pulse = 0.0003654 sec. 

required sampling time for 256 pulses = 0.0935424 sec. 

remaining time, for settling = 0.226 sec. 
25 decel sampling time = x 1.13 = 0.1057029 sec. 

remaining time, for settling = 0.214 sec. 

avg. finishes, if 40% decel, at 0.32 sec. 

avg. finishes, if no decel, at 0.308 sec. 

Ng: 



30 



max speed = 15,900 RPM = 265 RPS 
x 9.138 ppr = 2421.6 Hz 
x 2 = 4843.2 Hz 
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period per pulse = 0.0002064 sec. 
required sampling time for 256 pulses = 0.052884 sec. 
remaining time, for settling = 0.267 
decel sampling time = x 1.13 = 0.0597073 
5 remaining time, for settling = .260 sec. 



The resultant values stored in either the transmit buffer 1 65 
(Fig. 1) or the transmit buffer 205 may be utilized for a variety of 
purposes. One purpose may be to have the computer 1 58 provide a 

10 signal to an annunciator 219 (Fig. 1) to indicate to personnel 

operating the related rotary apparatus that excessive vibration has 
occurred, when that is the case. Similarly, the data may simply be 
recorded for consideration by maintenance personnel at a later date, 
such as after a plane has landed. 

1 5 In Fig. 1 , only a single vibration sensor 20 is shown. 

However, a second vibration sensor and corresponding charge 
amplifier, or even more of them, may be provided to the multiplexer 
22, and the apparatus of Figs. 1 and 6 will process all such vibration 
sensor outputs, one at a time. In the present embodiment, the 

20 apparatus may alternate between processing signals from the 

vibration sensor 20 and square wave test signals from the source 23. 
If additional vibration sensors are present, they may be processed 
contiguous with processing the vibration sensor 20, followed by 
processing of the square wave test signals, or the square wave test 

25 signals could, if desired, be interleaved with processing of a plurality 
of vibration sensors. Aspects of the invention, may, of course, be 
used with only one speed generator such as either Nf or Ng: and, if 
another channel of processing hardware and software is provided, 
the invention may be used with additional speed sensors. None of 

30 this is significant to the present invention. 
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It is intended that the Nf and Ng pulses determine the relative 
instantaneous angle of the fan and gas generator, and that they be 
utilized in a simple way to provide the sine and cosine of the angle of 
the fan and gas generator positions without any regard to an angular 
5 reference to a stationary position, or phase. 

It can be shown that the sine and cosine values for a fan 
speed (Nf) counter issuing 1 2 pulses per second, which is multiplied 
by four to avoid aliased errors due to sampling, repeats identically 
every 48 pulses. Therefore, only 48 sets of sine and cosine values 

10 for the 12 pulse per revolution fan speed sensor is required. This is 
accomplished in Fig. 1 by a modulo 48 counter 81 feeding a lookup 
table 82 of sines and cosines. The content of the lookup table 82 is 
partially shown in Fig. 5, where the sequence number is the output of 
the modulo 48 counter 81 . 

1 5 One example of an application for digital compensation of a 

non-integer number of pulses occurring with respect to a related 
event, such as a revolution of a rotating device, is in a case where 
the speed counter on the gas generator of a gas turbine engine 
provides 9.138 pulses per revolution, such as Ng herein. In this 

20 example, to avoid aliased errors due to sampling, the pulses per 

revolution (referred to as frequency, at times herein) is doubled, and 
the result is desired to be used for table lookup of sine and cosine 
values thereof. However, if these values are used, the values 
continuously change, out to one thousand times, before there is any 

25 cyclic recurrence which would be utilized in an integer number of 
pulses per revolution. Therefore, a lookup table would require over 
nine million entries. 

To handle the 9.1 38 pulses per revolution of the gas 
generator speed sensor, Ng, the invention first selects a more useful 
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non-integer number, in this case 9.125, and then determines that the 
difference in sine and cosine values between two times 9.138 and 
two times 9.1 25 are acceptably close, within a very small percent of 
each other, which in most utilizations is negligible and can be 
5 ignored. A small integer multiple, in this case four, of the number of 
pulses per revolution, in this case 18.25 yields a table having entries 
which repeat, in this case 73 entries (0-72). 

If there were 9 pulses per revolution, which are multiplied by 
2, to provide 1 8 pulses per revolution, the sine and cosine values 

10 thereof would repeat identically every 18 pulses, and also every 72 
pulses. According to the invention, assuming there are 9.125 pulses 
per revolution, and assuming those are doubled to 18.25 to avoid 
aliasing errors, then the sine and cosine values track almost perfectly 
and repeat after every 73 rd pulse. Thus, providing a table of sines 

1 5 and cosine values having 73 entries (that is, 73 values for successive 
input pulses) provides sine and cosine values which are adequately 
close in value to sine and cosine values that would be determined 
either by calculation or by an infinitely long storage device for the 
actual 9.138 pulses per revolution. 

20 Fig. 6 is fragmentary, showing a comparison at the first few 

counts, some middle counts and at the last few counts, illustrating 
that the values utilizing 18.25 definitely repeat after 73 pulses 
numbered 0 through 72 (the 73 rd pulse being zero) and also 
illustrating the degree of likeness. For instance, for the 65 th pulse: 

25 the value of the sine of 4 x 1 8.25 pulses per revolution is -0.38 while 
the value for the sine of 4 x 18.276 pulses per revolution is -0.35. 
The maximum error is on the order of 0.14%. 

Referring to Fig. 1 , the Nf signal from the fan speed generator 
43 is received on the line 44 and the frequency thereof is multiplied 
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by four in the frequency multiplier 45, as described with respect to 
Fig. 3. These signals are applied on the line 74 to the modulo 48 
counter 14 which simply counts the pulses and provides the count to 
the sine/cosine lookup table 82, The lookup table conforms to 
5 columns 1-3 of Fig. 4, referred to hereinbefore. The lookup table 82 
provides the exact sine of Nf on a trunk of lines 85 and the exact 
cosine of Nf on a trunk of lines 86. 

Similarly, the gas generator speed sensor pulses, Ng, are 
received on a line 89 from the gas generator speed generator 90, and 

10 multiplied by two in a multiplier 92, which may take the form of Fig. 

3, but only truncating by one bit, so n = 1 , 1/2 n is equal to one-half, 
and 2 n = 2. The doubled frequency signal is passed over a line 93 to 
a modulo 73 counter 95, the output of which is applied to a 
sine/cosine lookup table 96 which conforms to columns 1-3 of Fig. 5. 

1 5 The table 96 provides an approximate value of the sine of Ng on a 

trunk of lines 36 and an approximate value of the cosine of an Ng on 
a trunk of lines 37. 

In this case, the number selected as a substitute for the 
number of pulses per revolution is itself a non-integer, but it is one 

20 for which a relatively low harmonic is an integer. In the general case, 
the invention may be practiced utilizing substitute numbers which are 
either integers or non-integers. The invention may of course be 
utilized to assist generating other trigonometric values of signals 
indicative of rotation. 

25 In this case, a selected number, 9.125, is chosen to emulate 

input pulses having 9.138 pulses per revolution; the low harmonic of 
9.125 is the eighth harmonic, having a value of 73. Some other 
harmonic could be chosen (but in this case it would be wasteful of 
apparatus and processing time); other choices may be made, either to 
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get a more accurate approximation of the sine and cosine desired, or 
for other purposes. 

An exemplary, conventional three pole filter suitable for use 
as three pole filters 141, 146, 181, 186 in Figs. 1 and 4 is illustrated 
5 in Fig. 8. This filter will be performed by the state machine, for 

speed, by using constants which are derived in the computer 158 as 
a function of the period, T (as illustrated in Fig. 9), of the speed 
signals Nf, Ng on the lines 44, 89. The constants are illustrated in 
Fig. 9. 

10 When a three pole filter subroutine is called in step 141, 146, 

181, 1 86 of Fig. 4, the program will enter the three pole filter 
subroutine through a start point 220 and a first test 221 determines 
if this is the first pass or not since the start process 1 34 (Fig. 4). 
This will be determined by the cycle timing and synchronization 

1 5 provided by the computer 1 58 which will set a first pass flag in each 
filter. Initially, the first pass flag will be set in any of the filters, and 
a series of steps 224-226 will reset to zero those functions that are 
calculated and then used in subsequent passes. And then the first 
pass flag is reset in a step 228. A subroutine 229 will compute the 

20 quadratic filter output in accordance with the equation therein. Then 
a subroutine 230 will compute a lag filter output in accordance with 
the equation shown therein. The result of the filter is Y F (k), which is 
the filtered value of X(k), is transmitted to the state machine, Fig. 4, 
in a step 233. Then all of the factors which are utilized in 

25 subroutines 229, 230 are updated in a series of steps 235-237. 

Then, the program reverts to point 142, 147, 182, or 187 in the 
routine of Fig. 4 through a return point 240 in Fig. 8. 

The built-in test of the present invention is caused, by the 
computer 1 58, to be interleaved with processing vibration signals 
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from the vibration sensor 20. A series of different tests are 
performed, at different frequencies, so as to test the entire apparatus 
at frequencies that are pertinent to the fan speed and the gas 
generator speed. The constants of Fig. 9 are calculated separately 
5 for groups of built-in tests or individual built-in tests. The computer 
1 58 provides square wave signals at the desired sampling frequency. 
The amplitude of the square wave peak value (P) is determined for 
each of the individual test signals by the following relationship: 
P = 4/nA, where A is the amplitude of the desired sine wave peak. 

10 The multiplexer 22 therefore has a plurality of square wave test 
inputs, even though only a single input from the square wave test 
source 23 is shown. Because of the bandpass filtering in the signal 
conditioning circuits 26, and the three pole filtering described with 
respect to Fig. 8, the response to the harmonics of the square waves 

1 5 provided by the square wave test source 23 can be predicted for a 
properly operating filter (Figs. 2 and 8). Thus, one aspect of the 
present invention is testing of apparatus responding to signals having 
complex sinusoidal components by means of predicting the response 
of the fundamental and harmonic components of a related square 

20 wave. The results are indicative of the proper or improper 
functioning of the tracking filter. 

Thus, although the invention has been shown and described 
with respect to exemplary embodiments thereof, it should be 
understood by those skilled in the art that the foregoing and various 

25 other changes, omissions and additions may be made therein and 

thereto, without departing from the spirit and scope of the invention. 
I claim: 



